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We consider a general problem of laser pulse heating of spherical metal particles with the sizes 
ranging from nanometers to millimeters. We employ the exact Mie solutions of the diffraction 
problem and solve heat-transfer equations to determine the maximum temperature at the particle 
surface as a function of optical and thermometric parameters of the problem. The main attention 
is paid to the case when the thermometric conductivity of the particle is much larger than that of 
the environment, as it is in the case of metal particles in fluids. We show that in this case at any 
given finite duration of the laser pulse the maximum temperature rise as a function of the particle 
size reaches an absolute maximum at a certain finite size of the particle, and we suggest simple 
approximate analytical expressions for this dependence which covers the entire range of variations 
of the problem parameters and agree well with direct numerical simulations. 



PACS numbers: 44.05. +e, 42.62.Be, 82.50.-m 

INTRODUCTION 

The problem of laser pulse heating of absorbing par- 
ticles embedded in transparent liquid medium is im- 
portant for different applications, including stimulation 
of chemical reactions, laser sintering, selective killing 
of pathogenical bacteria or cancer cells, etc. A metal 
nanoparticle excited by laser pulse with the frequency 
close to its plasmon resonance efficiently converts electro- 
magnetic energy into thermal energy, with dramatic rais- 
ing of the temperature in the surrounding medium. As a 
result, this problem or its substantial parts were a sub- 
ject of many recent papers in physics, biology, medicine 
and chemistry (see, e.g., Refs. [THS] to cite a few). 

In spite of many experimental works on the subject, a 
theoretical understanding of this phenomenon is rather 
limited. In a standard approach, the problem is formu- 
lated as the study of heat transfer with a source (energy 
release in the particle) obtained as a solution of the cor- 



responding diffraction problem. Such a problem does not 
have simple exact analytical solutions, and different stud- 
ies employ either approximate analytical methods (valid 
for certain limiting cases only) or direct numerical cal- 
culations with the specified parameters. To the best of 
our knowledge, neither general solution of this problem 
applicable for a broad range of variations of the prob- 
lem parameters, nor suitable analytical expressions are 
available in the literature. 

On the other hand, for many applications it is highly 
desirable to obtain simple analytical expressions which 
describe the temperature at the surface of the particle 
in a broad range of the problem parameters, and for the 
particle sizes ranging from nanometers to millimeters. In 
the present paper we solve this problem. In particular, 
we obtain analytical expressions for the maximum tem- 
perature at the particle surface as a function of optical 
and thermometric parameters of the problem. Our re- 
sults cover a wide range of the particle sizes from a few 
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nanometers to millimeters, and they are confirmed by 
direct numerical solutions. 

The paper is organized as follows. In Sec. II we formu- 
late the problem and introduce four major spatial scales, 
namely the particle size, thickness of the skin layer, as 
well as the characteristic length of the heat diffusion in 
the particle and surrounding medium, which determine 
a variety of different cases discussed below. Our ana- 
lytical results and estimates are summarized in Sec. Ill 
for twelve different cases which include all possible com- 
bination of the four spatial scales. Section IV summa- 
rizes our numerical results which are found to be in a 
good agreement with the analytical predictions. Finally, 
Sec. V concludes the paper. 



PROBLEM FORMULATION 

We consider linear heating of metal particles by a single 
laser pulse with duration r. We assume that the particle 
has the thermometric conductivity Xp (typically Xp ~ 
0.1 — 1 cm^/s), and it is embedded into a fluid with the 
thermometric conductivity Xf (typically Xf ^ 10~^ — 
10~^ cm^/s), such that the condition Xp ^ Xf always 
holds, and the heat transfer is limited entirely by the 
heat diffusion. We neglect convection processes owing to 
large characteristic time required for the convection to 
arise and develop. 

Thus, the heating problem depends on the following set 
of parameters: the intensity of the laser pulse /, its dura- 
tion r, frequency w, optical and thermometric constants 
of the particle and environment. As for the intensity /, 
owing to the linearity of the problem, the heating rate 
is just proportional to /, so the results obtained for a 
given / may be easily recalculated for any other its value 
by a simple scale transformation, see below. As for the 
other parameters, if the particle is spherical with radius 
R and the thickness of the skin layer S (for metals in 
the optical region 6 ^ 10~^ cm), we can define just four 
characteristic spatial scales of the problem, namely R, 
S, yjxj^^ and ^JXp^- According to the general principle 
of the dimensional analysis, interplay between these four 
scales determines the entire variety of heating regimes of 
the particle. All these regimes are discussed below one 
by one. 



ANALYTICAL RESULTS 
Small particles 

For small particles, we assume that i? <C (5, so that the 
heating occurs in the entire bulk of the particle rather 
than on its surface. Interplay between other scales al- 
lows to consider several different cases which we analyze 
below. 



Case SI: i? <C (5 <C y/xp" ^ \^Xp^- this case 
the incident light penetrates into the entire particle and 
the energy release should be proportional to the particle 
volume, i.e. the absorption cross-section of the particle 
o'abs R^- Bearing in mind that the dimension of aats 
is cm^, it is convenient to write 



aabs = akR^, 



(1) 



where a is a dimensionless quantity, k stands for the 
wavenumber of the incident light in the environmental 
fluid: k = rifU) j c. Here is purely real refractive index 
of the fluid, w designates the frequency of the incident 
light and c stands for the speed of light in vacuum. 

In the simplest case of the Rayleigh scattering compar- 
ison of Eq. ([T]) with the well-known expression for Oabs 
for a spherical small particle |6j provides the following 
expression for a: 



(e' + 2)'+e": 



(2) 



where e' and e" stand for the real and imaginary parts of 
the complex relative dielectric permittivity of the particle 
(e = ep/n^). 

Next, we neglect spatial inhomogeneity of the heat 
sources, supposing that the volume density of the sources 
inside the particle is a constant equal to aabsI(t)/V : 
where I(t) is the power density of the laser pulse (in 
W/cm^) and V is the particle volume. The ground for 
this neglect is that the actual inhomogeneity in the heat 
sources is rather weak and it results even in weaker tem- 
perature inhomogeneities, owing to high rate of heat 
transfer in metals. It allows to replace the actual 3D heat 
transfer problem by its spherically-symmetric version. 

Moreover, in what follows we are interested in esti- 
mates of the maximal temperature of the particle surface, 
rather than in its exact calculations. For this reason we 
employ the spherically-symmetric problem formulation 
even for large particles, when the illuminated part of the 
particle obviously has temperature higher than that in 
shadow. 

Taking into account that inequality R <^ \/xf^ means 
the temperature field in the vicinity of the particle is 
quasi-steady (i.e., the term with the temporal derivative 
of T in the heat conduction problem may be neglected 
relative to the terms with the spatial derivatives) we ob- 
tain that within the framework of the approximations 
made the temperature field is described by the following 
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spherically-symmetric boundary- value problem: 



1 d 



'^P J.2 Qj. [ Qj. 



1 d f^^dT 



3al{t) 



47ri?3 
0, at r>R 



0, at r < i? 



dr \ dr 
T{t,R-0) = T{t,R + 0), 
dT\ fdT 
dr 



T — > at r — > oo, 



(3) 

(4) 
(5) 
(6) 

(7) 



where Kp and n / are thermal conductivity of the particle 
and environment, respectively and t in the dependence 
r(r, t) plays a role of a parameter. Here and in what 
follows T stands for the temperature rise from the room 
temperature, see Eq. ([tI. 

Integration of Eqs. |3|-([7| yields a parabolic temper- 
ature profile dX r < R and T = Ts{t)R/r at r > R with 
the surface temperature of the particle 



al{t) _ aR^kI{t) 



(8) 



Note, the temporal dependence of Tg coincides with that 
for I{t), so that the maximal temperature is achieved at 
the maximum of the laser pulse. 

Case S2: R ^ ^XfT ^ ^ <C ^Jx^- The case cor- 
responds to a quasi-steady field of the temperature both 
in the particle and in the fluid. From the point of view 
of heat conductivity, the case is identical to (SI), so that 
Eq. (|8| is vahd. 

Case S3: R <C y/Xf^ ^ \/Xf>^ ^ ^- The case is 
identical to (SI), and again Eq. (Is]) is valid. 



Large particles 

For large particle, we assume that i? ^ (5, so that the 
energy release occurs in a thin layer near the particle 
surface. 

Case LI (5 <C i? ^ \/Xj^ ^ \fX^- '^^^ ^^^^ 
absorption cross-section should be proportional to the 
square of the linear size of the particle, i.e.. 



chafes — T^R^Qab 



(9) 



where Qats stands for the dimensionless efficiency. 

Regarding the temperature field, owing to the condi- 
tion R ^ ^Xf' \/XpT it is still quasi-steady. The 
difference between the previous case is that now a good 
approximation is the energy release at the surface of the 
particle. Then, the temperature field inside the particle 
should satisfy the homogeneous Laplace equation. The 
only non-singular solution of this equation is a constant 
profile, so that T{r, t) at r < Ris reduced to T{t), where 
time t once again is regarded as a parameter. As for the 
temperature field in the environmental fluid, it keeps the 



same profile as that at R S, which eventually [bearing 
in mind Eq. (|9|] brings about the following expression 
for Ts-. 



alit) _ QabsRHt) 



AnRKf 



(10) 



Case L2: (5 <C ^Jxj^ <C i? <C ^/x^- In this case 
the field inside the particle is still quasi-steady, so that 
its r-dependence may be neglected, see above, case (LI). 
Regarding the field outside the particle, it is essentially 
t- and r-dependent. To determine we may employ 
the energy conservation law. For simplicity we consider 
a rectangular laser pulse with intensity /q and duration 
T. To a certain moment of time t < t the energy W 
absorbed by the particle is Gabslot- 

The absorbed energy is consumed to heat the particle 
to temperature Ts and to heat an adjacent layer of the 
fluid. The former requires the energy {A/3)TTR^CpPpTs, 
the latter AnR^2^/xJiCfPjTg/2, (to enhance accuracy of 
the estimate we have taken into account that the scale 
of a layer heated to time t is 2^/xjt [I] and replaced the 
proflle of the temperature in the heated layer by its mean 
value Ts/2). Here C and p stand for the specific heat and 
density of the particle (subscript p) and fluid (subscript 
/), respectively. 

Equalizing W to the consumed energy and considering 
the equality as an equation for unknown Ts, one easily 
derives 



Cabs lot 



lirR^CpPp + ATrR'^y 

Qabslot 

\RCpPp+A^^CfPf' 



(11) 



The obtained Tg{t) is a monotonic function of t, so the 
maximal temperature is achieved in the end of the laser 
pulse. Replacement t ^ t brings about the correspond- 
ing expression for the maximal temperature Ts(„jaj.). 

Case L3: (5 <C ^Jxf'' VXp^ ^ R- According to 
the employed problem formulation this case correspond 
to heating of inflnite compound space whose left semi- 
space has the thermometric properties of the particle, 
the right one those of the fluid, and energy is released 
at the boundary between the semi-spaces. This problem 
is exactly solvable [7]. The solution yields the following 
expression for the surface temperature: 



Tsit) 



(t/o 



VxpXft 



2R^7T^ Kpy/Xf + Kf^p 

ah VXpXft 

2TTy/n Kpy/Xf + >^-f^/Xp 



(12) 



Once again Tg occurs a monotonic function of time, so 
the maximal temperature is achieved in the end of the 
pulse, at t = T. 
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Case L4: ^Jxjr S <^ R >/Xp^- The surface 
absorption of light. The case is identical to (L2), so that 



Eq. (11) is valid. 

Case L5: yJxjT ^ ( 
identical to (L3), and Eq. (12) is valid 



The case is 



Other cases 

Case Ol: yjxs^ R <^ 5 y/Xp'T- From the view- 
point of the heat transfer the case is equivalent to (S2), 
but the absorption cross-section is described by Eq. ([T]). 

As a result 



Ts{t) 



Cabs lot 



^irR^a 



pPp + AnR^CfPfy/xfi 
akiQt 



(13) 



Case 02: 



6. The case is 



identical to (01), and Eq. (13) is valid. 

Case 03: ^X/t <C <^ R <^ S. According to 

the approximation used the energy is released inside the 
particle in a spatially-uniform manner. The case occurs 



analogous to(Ol), and again Eq. (13) is valid. 

Case 04: ^Jxjf <^ ^/Xp^^ ^ (5 ^ i?. In this case the 
energy release occurs in a surface layer with the volume 
AttR'^S. The condition ^Xp^ ^ ^ allows to neglect dis- 
tortion of the temperature inside the layer by heat con- 
ductivity. Thus, the released energy is consumed to heat 
the mentioned layer and the layer of fluid with thickness 
2y/xft. The energy balance yields the following expres- 
sion for Tc! 



Tsit) 



Cabs lot 



AnR^iCpPpS + CfPf^^) 

Qabslot 

A{CpPp5 + CfPfy/Yfiy 



(14) 



The maximal temperature is achieved in the end of the 
laser pulse, ai t = t. 

Thus the entire set of possible relation between the 
four length scales of the problem has been inspected. 



NUMERICAL VS. ANALYTICAL RESULTS 

Let us discuss the results obtained. Note that the de- 
pendence of the dimensionless quantities a and Qabs on R 
usually is very weak except for a narrow region centered 
about R ^ S, where the dependence aabs ^ R^ is replaced 
by (Tabs ^ R^ ■ To illustrate this point the dependencies 
Cabs{R) and the corresponding Qabs{R) for a spherical 
gold particle are presented in Fig. [T] The calculations 
are made according to the exact Mie solution [F and ac- 
tual optical constants of gold W. at the wavelength of the 
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FIG. 1. (color online) Absorption cross-section Gabs and the 
corresponding dimensionless efficiency Qabs for a spherical 
gold particle at A = 532 nm and (5(A) = 22 nm as functions 
of the particle radius R. Calculations according to the exact 
Mie solution. 
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FIG. 2. (color online) Comparison of the analytical and nu- 
merical results for the maximal temperature of a spherical 
gold particle with radius R in water heated by a rectangular 
laser pulse with wavelength 532 nm, lo = 2 ■ 10^ W/cm^ and 
T = 50 ns. [S <C ^yXfT, see cases (SI) and (L1)-(L3)]. 



incident light in a vacuum A — 532 nm. The correspond- 
ing value of 6 calculated as c/(cjn^„), where n'^^ stands 
for the imaginary part of the refractive index of gold at 
the given wavelength, is 22 nm. The cubic dependence 
CabsiR) [linear dependence of Qabs: i-e. independence of 
a of i?] at i? < (5 is seen straightforwardly. At R > 5 
the efficiency Qabs drops from 0.6 to 0.2 when R varies 
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FIG. 3. (color online) Same as in Fig. |2] at r = 50 ps. 
[\/x7t^ < 5 < -vAyr, see cases (S2), (01), and (L4)]. 
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FIG. 4. (color online) The same as that shown in Fig. [2] at 
r = 50 fs. [yx^ < S, see cases (S3), and (02)-(04)]. 



in three order of magnitude, i.e. the dependence QatsiR) 
is extremely weak, and our assumption aabs ^ does 
capture the main i?-dependence of Oabs in this area. 

It means that the main dependence of Tg(^jy^a^-^ on the 
particle size is given by the explicit i?-dependence of 



Eqs. (^, (10|-(14). In particular, for cases (i)-(iv) the 



quadratic growth of Tsf^^ax) with an increase in R is re- 
placed by the linear at i? ~ 5, then T, 



s{raax) 



reaches its 



maximum at i? \fXf^i declines with further increase 
in R and finally approaches to a constant at i? 3> JXv^- 
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FIG. 5. (color online) Density plot of the maximum of the 
surface temperature T^i^^^x) (arbitrary units) for a gold par- 
ticle in water irradiated by a rectangular laser pulse with 
Jo = 2 • 10^ W/cm^ as a function of the duration of the pulse 
T and the particle radius K. The dashed line corresponds 
to i? = 2,JxTf. Above this line the temperature becomes 
r-independent, see Eqs. (jsl, (10 1. 
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FIG. 6. (color online) The radius Rmax providing the absolute 
maximum of the surface temperature and the corresponding 
temperature Tg^Max) as functions of the laser pulse duration 
T for a gold particle in water irradiated by a rectangular laser 
pulse with lo — 2 ■ 10^ W/cm^ . The dashed line indicates de- 
pendence R — (a). The initial part of plot Tg(Max){'''), 
shown in panel (a), in Log-Log scale (b). 



For other case Tg(^„iax){R) is treated in a similar manner. 

Accuracy of the developed approach is illustrated by 
comparison of the obtained analytical expressions with 
numerical simulations of the corresponding spherically- 
symmetric version of the heat conduction equation for a 
gold particle in water presented in Fig. 2-4, where aabs 
is taken from the exact Mie solution for the gold particle. 
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The particle is heated by a rectangular laser pulse, whose 
intensity for definiteness is assigned the typical value Iq = 
2 ■ 10^ W/cm^, and various values of r (indicated in the 
figure captions). 

Note, that there are two competing mechanisms of 
maximization of the surface temperature. The first is re- 
lated to optics being associated with the local maximum 
of Qabs at i? ~ (5, see Fig. [T] The second is related to the 
heat transfer problem. It is associated with the change 
of the quasi-steady, i?-dependent temperature field in the 
vicinity of the particle to the essentially time-dependent 
temperature profile, which does not depend on R. The 
change occurs at i? « ^^fx}^- 

To understand the relative role of these mechanisms we 
plot the Ts{ynax) as a function of R and r, see Fig. [sj At 
every given R the Ts^max) increases monotonically with 
an increase in r until the latter reaches the values r ~ 
R? IXf. Then, Tg^jyjcix) becomes r- independent. 

Next, we perform the following calculations. For every 
given T we find such a value of i? = Rmax that maxi- 
mizes Ts(^jy^ax){R)i i-e., provides the absolute maximum 
of the surface temperature, which may be achieved for 
the given r. Then, Rmax and the corresponding temper- 
ature Ts(Max) = Ts(^rnax){Rmax) are plotted as functions 
of r, see Fig. [6j It is seen straightforwardly that for the 
problem in question heating for several degrees and high 
begins from r > 10^^ s, when the heat transfer mecha- 
nism prevails over the optic one. Let us stress that, as 
it has been already pointed out in Sec. II, owing to the 
linearity of the problem the temperature is just propor- 
tional to Jo, which allows easily recalculate the results 
obtained for a given value of Iq to any other its value. 

To illustrate how our results may be employed in vari- 
ous applications let us consider an important example of 
selective laser photo-thermal therapy of cancer. Thus, it 
was shown that 40 nm gold nanoparticle conjugated to 
certain antibodies and then incubated with both human 
oral cancer cells and nonmalignant skin cells were prefer- 
entially and specifically bound to the cancer cells. Next, 
the nanoparticle-labeled cells were exposed to a CW ar- 
gon ion laser at 514 nm. It was found that the malignant 
cells required less than half the laser energy to be killed 
as compared to the benign cells. The destruction of the 
cancer cells occurred owing to laser heating of the gold 
nanoparticles up to a certain threshold temperature [TU] . 

On the other hand, applying this approach in vivo to 
avoid unwanted effects one should minimize the expo- 
sure to the laser beam of the benign tissues and even the 
tumor itself, maximizing the rate of energy delivered to 
the nanoparticles and stored in the particles and their 
immediate vicinity. To this end a pulse laser should be 
employed. Application of our results indicates that re- 
duction of the pulse duration from infinity (for a CW 
laser) to a nunosecond scale (without change of /q) prac- 
tically does not affect the maximal temperature rise of 
the nanoparticles. Presumably, heating of the nanopar- 



ticles by such a pulse still should be fatal for the cancer 
cells bound to them, while the exposure of the rest of the 
tumor and the benign tissues will be reduced dramati- 
cally. 



CONCLUSION 

We have demonstrated that rather a complex problem 
of laser pulse heating of spherical metal particles embed- 
ded in a transparent fluid may be described by relatively 
simple analytical expressions, which provide the depen- 
dence of the maximum temperature at the particle sur- 
face Tg(^max) on the particle size and other parameters of 
the problem. We have demonstrated, by a direct com- 
parison with the numerical simulations, that these ex- 
pressions remain valid for any practically important size 
of the particle and duration of the laser pulse. More im- 
portantly, at fixed values of the material constants, the 
function Tg(^max){R) reaches local maxima at i? ~ (5 and 
R ~ \p(J^- For a gold particle in water considerable 
heating begins from r > 10^^ s, when the absolute max- 
imum of the temperature is achieved at i? w 2,JxfT. We 
believe these results not only give the most general so- 
lution of the problem but will also be useful to optimize 
the efficiency of laser pulse heating of nanoparticles in 
various problems of laser-matter interaction. 
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